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Abstract 

We confront the discovery of a boson decaying to two photons, as reported recently 
by ATLAS and CMS, with the corresponding predictions in the Minimal Supersym- 
metric Standard Model (MSSM) and the Next-to-Minimal Super symmetric Standard 
Model (NMSSM). We perform a scan over the relevant regions of parameter space in 
both models and evaluate the MSSM and NMSSM predictions for the dominant Higgs 
production channel and the photon photon decay channel. Taking into account the 
experimental constraints from previous direct searches, flavour physics, electroweak 
measurements as well as theoretical considerations, we find that a Higgs signal in the 
two photon channel with a rate equal to, or above, the SM prediction is viable over 
the full mass range 123 GeV < Mh < 127 GeV, both in the MSSM and the NMSSM. 
We find that besides the interpretation of a possible signal at about 125 GeV in terms 
of the lightest CT^-even Higgs boson, both the MSSM and the NMSSM permit also a 
viable interpretation where an observed state at about 125 GeV would correspond to 
the second-lightest CT^-even Higgs boson in the spectrum, which would be accompa- 
nied by another light Higgs with suppressed couplings to W and Z bosons. We find 
that a significant enhancement of the 77 rate, compatible with the signal strenghts 
observed by ATLAS and CMS, is possible in both the MSSM and the NMSSM, and 
we analyse in detail different mechanisms in the two models that can give rise to such 
an enhancement. We briefly discuss also our predictions in the two models for the 
production and subsequent decay into two photons of a CP-odd Higgs boson. 
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1 Introduction 



Two of the most important goals of the experiments at the Large Hadron Colhder (LHC) 
are to identify the origin of electroweak symmetry breaking, and to search for physics effects 
beyond the Standard Model (SM). The Higgs searches currently ongoing at the LHC (and 
previously carried out at the Tevatron) address both those goals. The spectacular discovery 
of a Higgs-like particle with a mass around Mh — 125 GeV, which has just been announced 
by ATLAS and CMS [T], marks a milestone of an effort that has been ongoing for almost 
half a century and opens a new era of particle physics. Both ATLAS and CMS reported a 
clear excess around ~ 125 GeV in the two photon channel as well as in the ZZ'^*'^ channel, 
whereas the analyses in other channels have a lower mass resolution and at present are less 
mature. The combined sensitivity in each of the experiments reaches about ~ 5cr. The 
observed rate in the 77 channel turns out to be considerably above the expectation for 
a SM Higgs both for ATLAS and CMS. While the statistical significance of this possible 
deviation from the SM prediction is not sufficient at present to draw a definite conclusion, 
if confirmed in the future it could be a first indication of a non-SM nature of the new state. 
The recent announcement follows the tantalising hints for an excess around Mh — 125 GeV 
that had already been reported both by ATLAS [2] and CMS [311]. The results of the SM 
Higgs searches at the Tevatron, based on the full dataset collected by CDF and DO have 
also just been announced [5], showing a broad excess in the region around Mh ~ 125 GeV 
that reaches a significance of nearly 3 a and would be compatible with a signal at about 
Mh ~ 125 GeV. 

The prime task now is clearly to study the properties of the discovered new particle and in 
particular to test whether the new particle is compatible with the Higgs boson of the SM or 
whether there are significant deviations from the SM predictions, which would point towards 
physics beyond the SM. The fact that the observed signal in the if — )■ 77 channel appears 
to be somewhat stronger than expected in the SM could be a first hint in this direction, 
although it is statistically not very significant up to now. This result nevertheless serves as 
a strong motivation for investigating possible alternatives to the SM where a possible signal 
in the 77 channel could be enhanced compared to the SM case. 

One of the leading candidates for physics beyond the SM is supersymmetry (SUSY), 
which doubles the particle degrees of freedom by predicting two scalar partners for all SM 
fermions, as well as fermionic partners to all bosons. The most widely studied SUSY frame- 
work is the Minimal Supersymmetric Standard Model (MSSM) which keeps the number 
of new fields and couplings to a minimum. The Higgs sector in particular contains two Higgs 
doublets, which in the CV conserving case leads to a physical spectrum consisting of two 
CP-even, one CP-odd and two charged Higgs bosons. 

Going beyond the MSSM, this model has a simple extension in the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM), see e.g. [7] for reviews. A particularly appealing 
motivation for considering the NMSSM is that it provides a solution for naturally associating 
an adequate scale to the n parameter appearing in the MSSM superpotential [8l[9]. In the 
NMSSM, the introduction of a new singlet superfield that only couples to the Higgs sector 
gives rise to an effective /i-term, generated in a similar way as the Yukawa mass terms of 
fermions through its vacuum expectation value. The new field must be a gauge singlet, 
since the parameter /i carries no SU{?>)c x SU{2)i x U{1)y quantum numbers. This effective 
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/x-term is linked dynamically to the electroweak scale. The additional degrees of freedom 
from the singlet add to the NMSSM particle spectrum. In the case where CV is conserved, 
which we assume throughout the paper, the states in the Higgs sector can now be classified 
as three CP-even Higgs bosons, hi {i = 1,2,3), two CP-odd Higgs bosons, aj (j = 1,2), and 
the charged Higgs boson pair H^. In addition, the SUSY partner of the singlet Higgs (called 
the singlino) extends the neutralino sector (to a total of five neutralinos) . 

The extended parameter space of the NMSSM as compared to the MSSM case gives rise 
to a rich and interesting phenomenology, in particular related to states in the Higgs sector 
that can be relatively light. As an example, Higgs-to-Higgs decays can be open (as in some 
of the benchmarks presented in [lO]), with potential consequences for Higgs searches at the 
LHC This non-trivial phenomenology of the NMSSM led in fact to the situation that it 
has not been possible to establish a general "no-loose" theorem for NMSSM Higgs discovery 
at the LHC [H] (similar conclusions were obtained for the case of the MSSM with complex 
parameters, see e.g. [I2].) Scenarios with a light CP-even Higgs boson with SM-like couplings 
to gauge bosons that however decays in an unusual way into a pair of very light CP-odd Higgs 
bosons ai (m^^ < 2mb), sometimes called "ideal Higgs scenarios" [13], have found particular 
attention. The case of a light CP-even Higgs in the NMSSM with suppressed couplings to 
gauge bosons has been discussed in [T3]. A summary of phenomenological results for the 
NMSSM Higgs sector can be found in [TSj. The presence of an additional neutralino can 
have implications for dark matter, which may in this case contain a sizeable singlet fraction 
(see [7] and references therein) or — like the new Higgs bosons — be very light |16] . 

Inspired by the spectacular discovery of a new particle in the LHC data that manifests 
itself in particular as a signal of at least (but possibly higher than) SM strength in the 77 
channel, we investigate the corresponding predictions in both the MSSM and the NMSSM, 
and compare them to the SM case. In particular, we evaluate the predictions for the produc- 
tion of a MSSM or NMSSM Higgs boson via gluon fusion, the main production channel at the 
LHC, followed by the decay into two photons, the channel with the largest significance in the 
results that have just been reported. As a main focus, we analyse potential enhancements 
of the production cross section times branching ratio over the corresponding SM prediction 
and we confront those predictions with the experimental data. We discuss in detail how an 
enhanced 77 rate can be realised in the MSSM, and which additional mechanisms for an 
enhancement can occur in the NMSSM. We find that sizable enhancements of the 77 rate 
are possible in both models, accompanied by or without a suppression of the WW^*^ decay 
mode. It is interesting to note that in both models the signal can be interpreted either as 
the lightest CP-even Higgs or as the second-lightest CP-even Higgs. The latter interpreta- 
tion would imply that there would be an additional Higgs boson present with a mass below 
125 GeV and suppressed couplings to gauge bosons. 

This paper is organised as follows: in Sect. |2]we give a short introduction to the NMSSM 
and specify our notation. In Sect. |3] we describe the more technical issues of the loop cal- 
culations involved: the framework for one-loop calculations in the NMSSM, the tools used 
for the numerical evaluation of observables, and we discuss constraints on the NMSSM pa- 
rameter space. The main results are given in Section |H We show the numerical results on 
loop-induced Higgs production and decays in the MSSM and NMSSM, which are compared 
to corresponding results in the SM (where relevant). Finally, the model predictions are con- 
fronted with the most recent results from the LHC Higgs searches, namely the discovery of 
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a Higgs-like particle, and the implications are discussed. Section |5] contains our conclusions. 



2 The NMSSM 

2.1 General considerations 

Like the MSSM, the NMSSM contains two Higgs doublet superfields Hi, H2, but in addition 
there is also a singlet scalar superfield 5*0 If one only considers the new interaction between 
the three scalar fields which gives rise to the effective //-term a Peccei-Quinn (PQ) symmetry 
is introduced. This may be avoided by considering a term cubic in the new singlet superfield, 
with a coupling strength k. This cubic term explicitly breaks the PQ symmetry, providing a 
mass to the scalar CP-odd (would be) axion. In order to reduce the interactions permitted in 
a general superpotential we assume a discrete Z3 symmetry. In this so-called Zs-symmetric 
version of the NMSSM the superpotential is assumed to be scale invariant (preventing terms 
linear or quadratic in S), thus it takes the general form 

^NMSSM ^ . ^^f)c ^ . jj^^c ^ ^ ^^^c ^ ^^^^ . ^ l ^^3^ 

where contraction of SU{2) indices is denoted by a dot product, i.e. Q ■ H = eijQ^H^ with 
eij = —eji and ei2 = 1. Assuming supersymmetry breaking by an unspecified mechanism, 
there are new soft breaking terms compared to the MSSM that need to be considered, 

^NMSSM ^ ^MSSM ^ ^2 |^|2 ^ ^^^^^^ . ^ 1 ^^^^3 + ^.c. (2) 

Ignoring the possibility of CV violation, as we do throughout this work, all parameters can 
be assumed real. 



2.2 The Higgs sector 

The explicit form of the Higgs potential is obtained from the Higgs F-terms, D-terms and 
the soft SUSY-breaking terms as Vh = Vp + Vd + Vgoit- The F-terms arise from the first 
derivative of the Za-invariant superpotential given in Eq. ([T]). As the singlet S is assumed 
to only couple to the Higgs, there is no coupling to the gauge bosons, and the D-terms are 
the same as in the MSSM. The soft breaking terms relevant to the Higgs sector are given in 
Eq. ([2]). Combining everything, we get the scalar Higgs potential of the NMSSM aqj 

Vh = ml\Hi\^ + mllH^l^ + ml\S\^ + \\eij{H2^Hlj) + kS^\^ + \X\^\S\' {\Hi\' + \H2\') 

+ (^\A^e,,{H2^H^,)S + ^-kA^S^ + c.c.^ + ^-{gl + gl) {\Hi\^ - \H2\^f + ]^gl\Hl^H2^? . 

(3) 

The indices {i, j} = {1,2} refer to the respective Higgs doublet component (summation over 
i and j is understood). The CP-invariant Higgs potential Vh contains the real soft breaking 

^Quantities with a hat denote superfields; fields without hats are their corresponding scalar components. 
^We use a notation that closely follows the one specified in [T7| for the MSSM. 
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parameters ml, m|, Ax, A^, as well as the superpotential trilinear couplings A and k, 
and the U{1)y and SU{2)i coupling constants gi and §2- 

After spontaneous symmetry breaking the two Higgs doublets can be expanded similarly 
to the MSSM. In components they are given by 

and the ratio of the two vacuum expectation values defines 

tan/3 = — . (5) 

The complex scalar singlet can likewise be expanded as 

S = Vs + ^{(ps + iXs) ■ (6) 

According to Eq. ([1]) a non-zero vacuum expectation value Vs for the singlet gives rise to an 
effective ix parameter (cf. the bilinear term W2 ~ ^iH2 ■ Hi in the MSSM superpotential) 

/icff = Xvs ■ (7) 

Since the scale for /lefr is set by Vs, a value at the electroweak scale is natural. Minimising 
the Higgs potential, the mass parameters m\ and m\ appearing in Eq. ^ are replaced by 
the vacuum expectation values Vi and V2 (or, equivalently, v = \/vf + ~ 174 GeV and 
tan/3). Similarly, the singlet soft mass m| can be expressed in terms of Vg- 
The bilinear part of the Higgs potential is then given by 




Vh = \ (01, 02, 0s) M^,^^ 02 +1 (Xi, X2, Xs) ^XXX X2 + (01 , 02 ) M^±^± ,+ + 




'2 



with the mass matrices M,^^,^, M^^^ and M(^±0± for the three CP-even neutral scalars 



(01,02,0s), three CP-odd neutral scalars (xi,X2,Xs), and two charged pairs (0]'',0^), re- 
spectively. At tree-level, the elements of the mass matrix M,^,^,^ for the CP-even Higgs 
bosons is given in the basis (0i, 02, 0s) by 

(M,^^,^) = M| cos^ /3 + 5/iefr tan /3, 
(M<^^<^)22 = M| sin^ /3 + 5/icfr cot /3, 

(M<^^<^)33 = cos /3 sin /3 + K^^^ {A^ + 4K/ieff) , 

/^off yd) 

(M000)^2 = (2^^f^ - Ml) cos /3 sin /3 - 5/ieff, 

(M000)^3 = \v [2/_icff cos /3 - (5 + -ft'/icff) sin /3] , 

(M000)23 = At; [2/_icff sin (3 - {B + Kfics) cos /3] , 

where we have introduced K = k/\ and B = Ax + Kfi^s, and Mz denotes the mass of the 
Z boson. 



4 



The CP-odd mass matrix is simplified significantly by introducing an effective doublet 
mass 

= 0^ (10) 

sm p cos p 

which corresponds to the mass of the single CP-odd Higgs boson in the MSSM limit. In the 



basis (xi, X2, Xs) this gives 




(Mxxx) 11 


= m\ sin^ (3, 


(Mxxx) 22 


= m\ cos^ /3, 


(Mxxx)33 


= {B + SK^es) COS 


(Mxxx) 12 


= m\ sin (3 cos /5, 


(Mxxx) 13 


= \v {B — 3Kfica) sin (3, 


(Mxxx) 23 


= \v {B — 3Kficf[) cos (3 



fill 



for the elements of the CP-odd mass matrix. 

The mass eigenstates in lowest order follow from unitary transformations of the original 
fields, 

= U^-(t], (a]=U^-(xX (c^±)=U^-('i]- (12) 




'A 




fai 


'2 


h 1 












'2 



The matrices ^ and t/*" transform the neutral CP-even, CP-odd and charged Higgs 
fields, respectively, such that the resulting mass matrices 

Mtt! = U"M,,,U"\ M^-J = U^M^^^U^' and Mf£^, = f/^M,.,.t/^^ (13) 

are diagonal in the basis of the transformed fields. The new fields correspond to the three 
neutral CP-even Higgs bosons hi, /12 and with rrih-^ < nih^ < rrih^, the two CP-odd Higgs 
bosons ai and 02 with < nia^, the charged pair and the Goldstone bosons G and G^. 
The tree-level mass of the physical charged Higgs boson is related to the CP-odd doublet 
mass through (with M^y denoting the mass of the W boson) 

Ml± =m\ + M^- \%\ (14) 

which shows that Mh± may be lower than its corresponding MSSM value. 



2.3 The neutralino and chargino sectors 

In addition to the extended scalar sector, the supersymmetric partner of the singlet Higgs 
field — the singlino S — leads to the presence of fifth neutralino in the NMSSM spectrum. As 
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a result, the neutralino mass matrix in the basis {B, H^, S) is given at tree-level by 



V 



]\/T 


n 
U 





M2 

























-9i 



92 



Vl 

V2 



9i 



-92 



V2 
V2 







V2 







-\V2 

-Xvi 



\ 



-Xvi 2K/iefr / 



It can be diagonalised by a single unitary (complex) matrix such that 

diag(m^o, m^o, . . . 



N*M^oN^ 



(15) 



(16) 



is real, positive definite, and with the mass eigenvalues ordered as rriy^o < m^o for i < j. 

Since no new charged degrees of freedom are introduced, the chargino sector of the 
NMSSM is identical to that in the MSSM. 



2.4 Squarks 

The squark sector of the NMSSM is unchanged with respect to the MSSM. Assuming minimal 
flavour violation and CV conservation, the mass matrix of the two squarks of the same flavour, 
ql and qr, is given by 

_ (Ml +ml + Ml cos 2/3(/| - Q,sl) m, X, \ 

^'"V ^'i^i + m2 + M| cos 2/3Q,4y'' ^''^ 



with 



Xq = Aq- /icff{cot /3, tan /?}, (1^ 



where {cot /3, tan/3} applies for up- and down-type squarks, respectively. In these equations, 
Mq^, Mq^ are (real) soft SUSY-breaking masses (as a consequence of SU{2) symmetry, Mq^ 
is equal for the members of an SU(2) doublet; in the numerical analysis below we assume a 
universal value Msusy = Mq^ = Mq^), while Aq is the soft SUSY-breaking trilinear coupling 
and yUeff the effective fi parameter as defined above. J| and Qq denote the weak isospin and 
the electric charge of the quark, respectively, and = \/l — with c„ = Mw/Mz- The 
squark mass eigenstates are obtained by the unitary transformation 

with 



Uq=( '^'% ''^'^ ), UqUl = l, (20) 

^ » — sm6'q COSb'q I ' 1 ' ^ ' 
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and the mass eigenvalues are given by 



+ M|^ + J|M| cos 2/3 



f211 



T J[Ml -Ml + Ml cos2/3(/| - 2Q,sl)Y + AmlXl 



3 Higher-order corrections in the (N)MSSM 
3.1 General considerations 

The investigations in the present paper are part of a bigger activity that we have undertaken 
in order to provide precise theoretical predictions for relevant observables in the NMSSM that 
can be used for comparing the NMSSM phenomenology with the SM, the MSSM and other 
scenarios of physics beyond the SM and for confronting the predictions with the available 
data. In order to investigate to what extent the Higgs search results can discriminate between 
different scenarios of physics at the Terascale, precise theoretical predictions both within the 
SM and possible alternatives of it are needed. In particular, if small deviations from the SM 
predictions are probed it is crucial to treat the considered model of new physics at the same 
level of precision to enable an accurate analysis and comparison. In the MSSM Higgs sector 
higher-order contributions are known to give numerically large effects (see, e.g., [TStITU]). and 
the same also holds for the NMSSM (see, e.g., [20] )• For many observables it is therefore 
necessary to include corrections beyond leading order in the perturbative expansion to obtain 
reliable results. It is planned that the results presented in this paper together with other 
results in the NMSSM Higgs sector [21] will be implemented into a new version of the Fortran 
code FeynHiggs [T71l22] - [2l] . which so far is restricted to predictions for Higgs physics in the 
SM and the MSSM. A few public codes already exist for numerical NMSSM calculations. 
The by far most widely used is NMSSMTools [25], which consists of subpackages to calculate 
the NMSSM spectrum, constraints on the parameter space, as well as Higgs decays and 
decay modes of sparticles [26]. Recently also an extension applicable to the NMSSM of the 
program SPheno [27] became available, which makes use of model implementations generated 
with SARAH [2S]. 

Since the NMSSM extends the MSSM in the Higgs and the neutralino sectors, differences 
to the MSSM are best probed in these two sectors. The processes playing the main role in 
the reported discovery at the LHC, production via gluon fusion and decay into two photons, 
are in fact processes that are particularly sensitive to possible deviations between the SM, 
the MSSM and the NMSSM. The one-loop predictions for those processes correspond to the 
leading-order contributions, which are IR- and also UV-finite without renormalisation (for 
a recent discussion of the renormalisation of the NMSSM Higgs sector, see [29]), so that the 
set-up mentioned above can immediately be applied for the investigation of the processes 
that are more important for NMSSM Higgs phenomenology at the LHC. After the initial 
announcements of an excess in the LHC Higgs searches by ATLAS [2] and CMS [3] there 
was already a considerable activity in the literature concerning possible interpretations of 
the results in various versions of the MSSM [30ti32] . and also in the NMSSM [33ti35] . We 
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comment on specific differences and similarities between those results and ours below. 

The calculation of loop diagrams, often involving a large number of fields, is a tedious and 
error-prone task if done by hand. This is true in particular for theories beyond the SM where 
the number of fields is significantly increased. For one-loop calculations, as will be the focus 
in the following, computer methods with a high degree of automation have been devised to 
simplify the work. However, most of the available tools so far have focused on calculations 
either in the SM or the MSSM. In order to facilitate loop calculations in the NMSSM, it 
is useful to employ the well-established public tools FeynArts ^36|, FormCalc [371|38] and 
LoopTools [371[39], with which one- loop calculations can be carried out with a high degree 
of automation. The program FeynArts [36] can be used to generate and draw the Feynman 
diagrams to a given order for the process under study, based on the information about the 
particle content and interactions that is supplied in a so-called model file. From the Feynman 
rules a mathematical expression for the corresponding amplitudes is generated. For one-loop 
amplitudes, the analytic simplifications, trace evaluation, tensor decomposition, etc. can then 
be carried out by FormCalc [371|38], which combines the speed of FORM |l0] with the more 
user-friendly interface of Mathematica. If the necessary information about the particles 
involved in the considered process (mass eigenvalues, mixing angles, etc.) is provided, the 
numerical evaluation of the FormCalc output can be carried out with LoopTools [37tl39]. 
In the standard distribution, FeynArts contains model files for the SM (including several 
variants), a general two-Higgs-doublet model, and an implementation of the MSSM. 

With the goal of treating the NMSSM at the same level of accuracy as the MSSM, we are 
currently developing a framework that enables the computation of one-loop processes in the 
NMSSM in a highly automated way. This consists on the one hand of the input on the particle 
content, the interaction vertices, etc. that is needed for use with the packages FeynArts 
and LoopTools, and on the other hand of an appropriate renormalisation prescription for 
the NMSSM. In the present paper we will focus on loop-induced processes, for which no 
renormalisation is needed. A detailed discussion of our renormalisation prescription for the 
NMSSM will be presented elsewhere. 

3.2 Implementation 

As a first step we have compiled a new FeynArts model file for the NMSSM. Since it is fore- 
seen that the same framework will be applied for several calculations beyond those presented 
in this paper, we describe our NMSSM implementation in some detail here. The basis for the 
model file itself — defining the particle content and interactions of the NMSSM (as described 
in Sect. |2]) in a general gauge — is generated with the help of the program SARAH [28] . 
This program can be used to generate FeynArts model files, as well as output for many other 
programs, for any supersymmetric theory starting from its superpotential. For consistency 
checks, we also use an independent NMSSM model file generated with FeynRules [H]. Start- 
ing from the output of SARAH we have introduced the standard nomenclature of FormCalc 
to activate its internal MSSM simplifications, applied unitarity relations to mixing matrices 
and couplings, and implemented some further improvements. These modifications, besides 
greatly improving the speed at which FormCalc performs one-loop calculations of NMSSM 
amplitudes, are essential for instance for verifying the cancellation of UV divergences at the 
algebraic level. 
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To enable the numerical evaluation of observables, such as decay widths and cross sec- 
tions, the analytic amplitudes of FormCalc can be exported to Fortran code. These must be 
supplemented with a so-called driver program to compute the necessary quantities (masses, 
mixings, etc.) from the fundamental parameters of the theory. The driver codes also provide 
standard facilities for numerical integration and the evaluation of master one-loop integrals 
through LoopTools. We have developed such a driver program for the NMSSM, which in 
its present state allows for Higgs and sparticle masses to be calculated either at tree-level 
(following Sect. [2]), or using NMSSMTools (here we used version 2.3.5) [25] linked through a 
custom interface. 

The NMSSM driver also offers the possibility to impose restrictions on the NMSSM pa- 
rameter space resulting from the evaluation of various experimental or theoretical constraints. 
For instance, the constraints implemented in NMSSMTools can be accessed, and direct con- 
straints on the extended Higgs sector are available through an interface to HiggsBounds ^2] . 
More details on the different constraints and how they are evaluated is given in Sect. 13.41 
below. 

3.3 Verification 

We have performed several tests on the model file to verify the NMSSM implementation, in 
particular for the Higgs sector which is most relevant for the present work. The analytical 
expressions for the Feynman rules for the interaction vertices of the NMSSM obtained from 
SARAH have been compared to the independent FeynRules output. They have also been 
compared (analytically) in the MSSM limit to the corresponding vertices in the default 
MSSM implementation distributed with FeynArts. 

A number of tree-level processes have been analysed numerically — including the decays 
of Higgs bosons and neutralinos — to test the mixing properties of the singlet state in the 
NMSSM. Comparing these to the resuhs of NMSSMTools and NMSDECAY [26], we find overall 
good agreement with those previously obtained results after correcting for differences due to 
QCD corrections and the running of gauge couplings. 

A further, extensive, and non-trivial test of the working NMSSM implementation is 
provided by the results for the processes that are induced at the one-loop level in the NMSSM. 
We have evaluated C(50) 1—7-2 processes and (9(100) 2—7-2 processes of this type and 
checked them successfully for their UV- and IR-finiteness. This includes the Higgs production 
and decay modes that are phenomenologically most relevant at the LHC, which are described 
and analysed numerically in detail below. As an example we show in Fig. [1] generic Feynman 
diagrams contributing to gg — >■ hi (upper row), where hi denotes any neutral CP-even Higgs 
boson in the (N)MSSM, as well as to hi — ?- 77 (lower row). These types of diagrams have 
been evaluated and used for our numerical analysis described in Sect. HI 

3.4 Constraints on the parameter space 

Before moving on to our numerical analysis, we briefly discuss the various phenomenologi- 
cal constraints which exist on the parameter space of the NMSSM. As already mentioned, 
our computational framework allows for all these constraints to be evaluated and applied 
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Figure 1: Generic diagrams contributing to gg — )• hi (upper row) and to hi — )■ 77 (lower 
row), where hi denotes any neutral CP-even Higgs boson in the (N)MSSM. 



in connection with the calculation of arbitrary NMSSM observables. In this way reliable 
predictions can be obtained which are in agreement with the present experimental results. 

Theoretical constraints 

Constraints originate from the requirement of a viable physical minimum of the Higgs po- 
tential. The physical minimum, with non-vanishing vacuum expectation values (vevs) for 
the two Higgs doublets Hi and H21 should be lower than the local minimum with vanishing 
vevs for the two Higgs doublets. Furthermore the physical vacuum should have a non-zero 
singlet vev to be able to generate the parameter. In the NMSSM, the correct pat- 
tern of a symmetry breaking absolute minimum is ensured approximately by the condition 
A\ > 9m| [13]. For each point considered in the NMSSM parameter space, we verified 
numerically (using NMSSMTools) that the Higgs potential is bounded from below and stable. 

Another requirement is that there be no Landau pole for any of the running couplings 
A, K, yt and yb below Mgut ~ 2 x 10^^ GeV. The renormalisation group equations for the 
NMSSM are known to two-loop order [H]. The constraint of perturbativity up to a very 
high scale restricts the range of A and n. Values of these parameters in the perturbative 
regime at the GUT scale lead to comparably small values at the weak scale, which may be 
combined to give the approximate upper bound [9] 

A^ + < 0.5. (22) 

In the parameter scan below we choose the ranges for A and k to respect this limit. 
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Bounds from direct Higgs searches 

The hmits from Higgs searches at LEP, the Tevatron, and the 2011 LHC data put bounds 
on the Higgs masses and couplings. The LHC bounds have now been updated [1] and allow 
a SM-like Higgs boson only in a small window around ~ 125 GeV. As strict bounds on the 
parameter space we will consider the data from LEP and the Tevatron as well as the LHC 
data presented in 2011 (referred to as LHC2011) and comment on the most recent limits as 
presented in pQ separately. 

In the NMSSM, or any other theory with a Higgs sector different from that in the SM, the 
lower limit on the Higgs boson mass of Mh > 114.4 GeV [15] from the LEP Higgs searches 
does not apply generically. The same is true for the limits from the searches for a SM-like 
Higgs at the LHC2011, which essentially rule out the mass range 127 GeV ^ ^ 600 GeV 
for a SM-like Higgs. As has been demonstrated in various benchmark scenarios in the 
MSSM [12] and the NMSSM [17], allowed parameter regions for Higgs bosons of those ex- 
tended Higgs sectors exist both below the LEP limit and above 127 GeV. In order to test 
whether a given point in the (N)MSSM parameter space is allowed or ruled out by the Higgs 
searches at LEP, the Tevatron, and the LHC2011 one therefore needs to confront the predic- 
tions of the model with the available cross section limits in the various search channels at 
each collider. For this purpose we make use of the code HiggsBounds [12] (version 3.6.1- 
beta). In order to obtain the correct statistical interpretation in terms of a 95% CL exclusion 
limit, HiggsBounds uses the input provided for the model under consideration (in the case of 
our analysis effective couplings and partial widths of the Higgs bosons of the MSSM and the 
NMSSM) to determine for each considered parameter point the channel that has the highest 
expected sensitivity for an exclusion. Only for this particular channel the theory prediction 
is then compared to the observed experimental limit, which determines in a statistically 
consistent manner whether the parameter point is allowed or excluded at the 95% CL. The 
version of HiggsBounds used for our analysis includes in particular limits from CMS on an 
MSSM Higgs decay into a tau pair [IS], which have a strong impact on the allowed region in 
the (Ma, tan /3) plane! Results from /i — >■ 77 [SUlCTj and h — )■ WW channels [5^I5B] (among 
others) are also implemented]!] Due to their importance to our analysis, we will also display 
these limits separately below. 

Bounds from direct searches for SUSY particles 

Constraints also exist from direct searches for supersymmetric particles at LEP, the Tevatron, 
and the LHC. The least model-dependent limits are the ones from LEP. In particular the 
limit on the lightest chargino mass, m^± > 94 GeV [SI], applies to both the MSSM and the 
NMSSM and restricts the parameter /i (/ieg) of the (N)MSSM to values above about 100 GeV. 

■^Using the latest limits from the CMS update |49] in this channel would lead to only minor changes in 
the results of our analysis. 

""^Strictly speaking, these results can only be directly interpreted in models with SM-like Higgs production. 
Even if this does not hold exactly, we have verified for the scenarios considered below that the ratio between 
the contributions from gluon fusion to vector boson fusion do not deviate from the SM case to the degree that 
vector boson fusion becomes the dominant production mode. Since we arc mainly interested in investigating 
scenarios with an enhanced rate, we take a conservative approach and apply these bounds for maximal 
exclusion. 
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For the squarks of the first two generations the LHC mass limits are most stringent [55] , 
but have a certain model dependence. The ATLAS and CMS collaborations set limits above 
1 TeV in particular "simplified models" for the masses of first and second generation squarks 
and for the gluino; slightly weaker bounds are derived in the constrained MSSM. For the 
third generation squarks, which are most important for the radiative corrections in the Higgs 
sector, the bounds are much weaker and more model dependent. This means in particular 
that the presence of large stop mixing (as favoured by a relatively high value of the Higgs 
mass of about 125 GeV, see for instance Ref. [30]) is not excluded at present. 

The anomalous magnetic moment of the muon 

A significant deviation exists between the experimentally measured value of the anomalous 
magnetic moment of the muon, a'^^^ and the theoretical prediction in the SM. Numerically 
this deviation amounts to Aa^ = a^^^^ — a^ = (30.2 ±8.8) x 10~^° [56|[57]. which corresponds 
to more than a 3 a effect. Employing new physics contributions to account for this deviation 
leads to bounds on the model parameters [58|[59]. The dominant contributions to in the 
NMSSM are known including leading corrections up to the two- loop order [59]. For the 
numerical evaluation we use NMSSMTools. As the 2cr allowed range for the NMSSM-specific 
contributions Aaj^^ssM _ ^nmssm_^sm ^ ^l x 10"^ < Aa^^ssM < 4 >< ^Yiich 

includes a theory uncertainty on the SUSY evaluation corresponding to 2.0 x lO"^'' added in 
quadrature to the uncertainty quoted above. We note that, similarly to the /i parameter in 
the MSSM, a positive value for fi^s is strongly favoured when is included as a constraint. 

Flavour physics 

A recent analysis and summary of flavour physics constraints on the NMSSM parameter 
space has been presented in [1?]. In the present setup we use NMSSMTools (version 2.3.5) to 
evaluate the NMSSM theory predictions. The corresponding experimental limits are listed 
in Table [H Parameter-dependent theory uncertainties are added linearly to the intervals 
shown in the table before evaluating exclusion. 

In theories with minimal flavour violation (MFV), which we are investigating here, the 
strongest constraints from flavour physics can usually be derived from S-physics observables 
such as BR(i? — > Xs'y), BR{Bs — )■ fj.~^fi~), BR(i?„ — )■ r"*"//,-), or from the mass mixings AMg, 
AMd [EniEl]- In addition to the MSSM-type contributions from sparticles and non SM-like 
Higgs states (charged and neutral), the NMSSM may also be further constrained in the case 
where a very light CP-odd Higgs boson is present [62]. 

It has been argued that NMSSM specific contributions, i.e. contributions that go beyond 
the MSSM, to BR(_B — )■ Xs'j) are negligible [61]. Besides the charged Higgs boson contri- 
bution, which typically yields large effects for a relatively small Mh±, we therefore include 
also the relevant loops of charginos, gluinos, as well as (non-singlet) neutralinos. Depending 
on the SUSY parameters, large cancellations between different contributions are possible, 
and any of the contributions mentioned above can become dominant. If a specific bench- 
mark scenario is chosen for the SUSY parameters in the other sectors, it cannot be expected 
that such a cancellation takes place for all values of the Higgs sector parameters. However, 
a comparably small shift in the underlying SUSY scenario would typically be sufficient to 
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Observable 


Exp. lower limit 


Exp. upper limit 




BR{B X,7)e^>i.6 GeV 


3.03 X 10-^ 


4.07 X 10-4 


[63] 






1.1 X 10-^ 


m 




0.79 X 10"^ 


2.45 X 10-4 


m 




17.53 ps-i 


18.01 ps"^ 


[M] 




0.499 ps-i 


0.515 ps-i 


[63] 



Table 1: Experimentally allowed ranges at the 2 a level used for the flavour physics observ- 
ables. 



give rise to large compensations between the different contributions, which could bring the 
prediction into agreement with the experimental result. 

Both the processes BR(_Bs — )■ and BR{Bu — t+z/,-) are also very sensitive to 

new physics effects, and can mainly exclude regions of parameter space with large tan/3. 
But while unacceptably large contributions to BR{Bs — )■ fJ''^lJ'~) can be avoided in a similar 
fashion as for BR{B Xg'^) by a slight change of the other parameters of the SUSY 
scenario, such a compensation in general does not work for BR{Bu t^Vt). This is due to 
the fact that this process involves charged Higgs exchange at tree-level. As a consequence, 
this process provides a similar exclusion power in the (M/^±, tan (5) plane for the MSSM and 
the NMSSM. 



4 Numerical analysis 

In this section we analyse numerically the phenomenologically important loop-induced Higgs 
decays of the neutral CP-even Higgs bosons to two photons, 

/i, ^77 (^ = 1,2,3) . (23) 

For completeness, we also discuss the decays of the CP-odd Higgs bosons, 

«i^77(j = l,2). (24) 

We investigate in particular, taking into account the existing constraints on the parameter 
space discussed above, to what extent the phenomenology of Higgs decays into two photons 
can differ in the MSSM and the NMSSM from the SM case. We choose a baseline MSSM 
benchmark scenario, fixing the soft SUSY-breaking parameters according to: 

MsusY = Mg-, = M,-^ = 1000 GeV, 

= Mi^ = 250 GeV (to comply with Aa^), 

At = Ab = Ar, 



5 1 

i\ /I 

2, 



Ml = ^tan^^wMz ^ ^M; 
Ms = 400 GeV 

m-g = 1200 GeV. (25) 
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While in the MSSM the tree-level Higgs sector can be specified by the two parameters Mh± 
(or Ma) and tan/3, the NMSSM Higgs sector has larger freedom and requires additional in- 
put. We choose the following set of parameters to describe a point in the NMSSM parameter 
space: 

Mh±, tan/3, \ K = k/X, . (26) 

The parameter Mh± here in principle plays the same role as in the MSSM. However, since 
we employ NMSSMTools to calculate the Higgs masses, the input Mh± is not defined in the 
on-shell renormalisation scheme, and must be understood as a tree-level input mass which 
can be directly translated to a value for Ax using Eqs. ( fTOl) . (fT4l) . The calculated physical 
Mh± (including the higher order corrections) will therefore in general not be identical to 
the input valuelj The parameter /Xefr, which we take as input, is in close correspondence to 
the n parameter in the MSSM. Together with A it determines the value of the singlet vev 
through Eq. ([7]). 

Below we are going to perform scans over the NMSSM parameter space. For the com- 
parison between the NMSSM, the MSSM, and the SM we classify the resulting scenarios as 
follows: 

(i) General NMSSM 

This refers to any scenario in the NMSSM specified by the full set of parameters given 
by Eq. (126!) and the MSSM parameters in Eq. (125|) (where we also defined our default 
settings and values). 

(ii) The MSSM limit 

To recover the MSSM, we start from an arbitrary NMSSM scenario and take 

A -> 0, K -> 0, K ^ k/\ = constant. (27) 

All other parameters (including fics) are held fixed, which corresponds to Vg — ?■ oo. 
In this limit the Higgs doublet (and neutralino) couplings become MSSM-like, and 
the couplings of the singlets vanish. It should be noted, however, that the decoupled 
singlets do not necessarily correspond to the heaviest states hs, 02, and X5 even in this 
limit. In the MSSM limit the parameters Mh± and tan/3 correspond exactly to their 
MSSM counterparts. 

(iii) Decoupling (SM) limit 

On top of the condition given by Eq. fl27|) for the MSSM limit, we define decoupling 
to the SM case by taking the additional limit 

Mh± > Mz. (28) 

This leads to heavy, nearly mass-degenerate, doublet Higgs states and leaves one light 
Higgs boson, H, with SM-like couplings to vector bosons and fermions. The singlet 
states are not affected by this limit; they remain decoupled with finite masses (which 
may or may not be lower than Mh)- 

^ This feature would be avoided with an on-shell renormalisation of M^±, see e.g. [T71I29] . 
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Parameter 


Minimum 


Maximum 


^SUSY 


750 


1500 


M2 ~ 2Mi 


200 


500 


At = Ah = Ar 


-2400 


2400 


^^ 


200 


3000 


Ma 


100 


600 


tan /3 


1 


60 



Table 2: Parameter ranges for the MSSM scan. All parameters with mass dimension are 
given in GeV. 

(iv) SM+singlet limit 

When the doublet decoupling condition Mh± ^ Mz is fulfilled for points with finite 
non-zero A, n (i.e., values that differ from the MSSM limit) we speak of the SM+singlet 
limit. This name is appropriate, since the low mass Higgs spectrum in this scenario 
consists of two CT^-even and one CP-odd degree of freedom. 

4.1 Decays of CP- even Higgs bosons in the MSSM 

Before we proceed to the NMSSM case, we study the two photon decays of the two CP-even 
Higgs bosons, h and if, in the MSSM and compare to the SM. For the numerical evaluation 
in the MSSM we use the code FeynHiggs (version 2.8.6) [T7 tl22] - [2^ . which is also used to 
evaluate SM quantities given in this section. To study interesting regions of the MSSM 
parameter space, where differences in the di-photon channel between the MSSM and the 
SM can occur, we perform a random scan over the parameter ranges given in Tab. |2j The 
remaining MSSM parameters are kept at their 'benchmark' values specified in Eq. fl25|l . 

It should be noted that we allow for comparably high values for /i; this is relevant for the 
possible size of some of the effects that we will discuss in detail below. However, such large 
values of /i, together with large values of tan /3, can lead to parameter combinations that show 
a non-stable behaviour in perturbation theory. In order to avoid parameter combinations 
that result in unacceptably large two-loop corrections in the evaluation of the Higgs boson 
self-energies and related quantities, we implement an upper limit on the corrections to the 
elements of the Z matrix (see Ref. [IZ]). Comparing the one- and two-loop values of the 
respective diagonal elements, we require the following condition for the light CP-even Higgs, 
I |2-2-ioop| _ i^'J^-ioopj |y'|2'^-i°°P| < 0.25, and analogously for the heavy CP-even Higgs with the 
replacement Zn — Z22. We found that this upper bound is effective for avoiding parameter 
regions that are deemed unstable under higher-order corrections. 

In Fig. [2] we show V{h — 77) in the top left and BR(/z — )■ 77) in the bottom left plot 
as a function of M^. The corresponding plots for if — )■ 77 are given in the right column. 
The colour coding is as follows: all points in the scan which are allowed by the theoretical 
constraints and the direct search limits for sparticles [51], as discussed above, are plotted in 
grey. Points which are also allowed by direct Higgs search limits (from HiggsBounds 3.6.1, 
i.e. including LHC2011) are shown in blue (on top of the grey points). Finally, points which 
fulfil additionally the constraint from — and BR(fe — )■ 37) (both are here calculated with 
FeynHiggs) are plotted in black. The red (solid) curve in Fig. [2] shows the corresponding 
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Figure 2: Results from the MSSM parameter scan for the partial widths r(/i, H — 77) of 
h (left) and H (right), and the corresponding branching ratios. The full result of the scan 
(all points allowed by the theoretical constraints and the direct search limits for sparticles) 
is shown in grey. The blue points are compatible with the direct Higgs search limits (from 
HiggsBounds 3.6.1, i.e. including LHC2011), while the black points in addition give a result 
in agreement with {g — 2)^ and BR(6 — )■ 57). The solid (red) curve shows the respective 
quantities evaluated in the SM. 



SM result with Mhq^ set equal to the corresponding MSSM Higgs mass. It should be noted 
that here (and in all the following plots) different densities of points appearing in different 
regions have no physical meaning, as the point density is related to the specific procedure 
chosen for the sampling of the SUSY parameter space. 

We first focus on the light CP-even Higgs boson, h, decaying into two photons. The 
extra particles in the MSSM yield additional loop contributions, which can both lower and 
raise r(/i — )■ 77) compared to the SM case. For Mh < 114.4 GeVl most of the scenarios 
where T{h — > 77) > T{Hsm ll) are ruled out by the direct Higgs search hmits, but we 
also find allowed points in this region. For those h couples with about SM strength to gauge 

^We neglect here, and in the foUowing, the theory uncertainty of the Higgs boson mass evaluation, which 
for the light Higgs boson should be roughly at the level of 2 — 3 GeV [24] . 
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bosons, but is nevertheless not excluded due to a (much) suppressed coupling to b quarks, 
which weakens the corresponding LEP limit. The fact that an enhanced rate for the decay 
of a Higgs boson to two photons is possible even below the LEP limit provides a motivation 
to extend the LHC Higgs searches to this region. In the following we focus on the mass 
region above the LEP limit. There we find scenarios in which r{h — )■ 77) is enhanced by 
up to ~ 70% with respect to the SM. On the other hand, as can be seen from the lower left 
plot in Fig. [21 the BIi{h — )■ 77) can be enhanced by a factor ~ 3 over the SM in the same 
mass range (due to a suppression of the bb decay mode as discussed in more detail below). 
It is interesting to note that for the points that are allowed by all constraints the maximum 
enhancement of the branching ratio occurs around ~ 125 GeV. One can finally observe 
that no SM values are reached for Mh ^ 114.4 GeV, reflecting the fact that a SM-like Higgs 
boson is ruled out by the LEP Higgs searches. 

The corresponding results for the heavy CP-even MSSM Higgs boson are shown in the 
right column of Fig. [2J For Mh ^ 130 GeV we find viable points with a BR slightly larger 
than for a SM Higgs boson. For larger values of Mh one can see the behaviour expected 
from the decoupling properties of the MSSM, i.e. T{H ^ 'j'j) and BR{H — 77) are both 
suppressed with respect to the SM, with the level of suppression increasing with Mh- 

In order to investigate the phenomenology at the LHC, besides the branching ratio of 
course also the Higgs production cross section — which in general will also be modified in a 
model of physics beyond the SM — has to be taken into account. The combined enhancement 
or suppression over the SM for a process pp ^ hi ^ X can therefore be summarised in the 
ratio 



. ^ a{pp hi) X BRjhj X) 
^ a{pp ^ Hsu) X BR(/fsM ^ X) ' 



(29) 



In the MSSM hi denotes either h or H, while in the NMSSM (discussed in the next subsec- 
tion) hi {i = 1 . . .3) can be any of the three CP-even Higgs states. If the Higgs production 
cross section is dominated by a single mechanism, such as gluon fusion which is often the case 
at the LHCfl, a common approximation is to use instead of a{pp — )■ hi) the parton-level cross 
section a{gg — > hi). Neglecting the differences in kinematics, the decay width T{hi — > gg) 
has the same dependence as algg hi) on the couplings of the involved particles, and the 
dominant higher-order QCD corrections are expected to cancel out in the ratioj^ Making use 
of this approximation, Eq. f l2^ can be expressed as 

r.hi ^ r(/i, ^ gg) X BRjhj ^ X) ^ T{h, ^ gg) x Tjhi ^ X) x T^otiHsu) , . 

T{HsM-^gg)y<BR{HsM-^x) t{Hsm^ gg) y<r{HsM^ x) xTUh^y ^ ' 

This definition will be used to calculate R^i^ and in the MSSM (and also in the NMSSM 
below) . 



^ We have checked that for the relevant regions of parameter space discussed below the gluon fusion 
production cross section always strongly dominates over the associated Higgs boson production from bottom 
quarks. 

^ Non-negligible differences are mainly expected if the bottom loop contribution to hi — > gg dominates 
over the top loop contribution. In the case of the light CT'-even Higgs boson can happen for very low 
Ma and moderate to large tan/3 values, whereas in the case of the heavy CP-even Higgs boson this can 
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Mh (GeV) Mh (GeV) 

Figure 3: Results from the MSSM parameter scan on the ratios R^^ for the hght CP-even 
Higgs boson h (left column) and R^^ for the heavy CP-even Higgs boson H (right column). 
The plots are displayed both on a logarithmic scale (upper row) and on a linear scale (lower 
row). The colour coding for the scan points is the same as in Fig. |2l The green lines are the 
corresponding limits from data presented in 2011 data from ATLAS [50] (solid) and from 
CMS [51] (dashed). The red lines are the new limits from ATLAS (solid) and CMS (dashed) 
taken from [1]. 

The results for R^^ and R^^ are shown in Fig. [3l with the same colour coding as in Fig. |2l 
In order to make the results better visible we display them twice on a logarithmic scale (upper 
row) and on a linear scale (lower row). The green curves in Fig. [3] show exclusion limits in 
the di-photon channel at 95% CL from data presented in 2011 data from ATLAS [50] (solid) 
and CMS [51] (dashed) @ The red lines are the new limits from ATLAS (solid) and CMS 
(dashed) taken from The exclusion limits from ATLAS and CMS are displayed here 

happen for larger Ma and tan/3 ^ 5. Our results therefore exhibit an additional uncertainly in this part of 
the parameter space. Additional loop contributions from SUSY particles, while taken into account in our 
calculation, are usually subdominant and of lesser importance in this context. 

^The 2011 exclusion limit from CMS in this channel was updated in [67 . Including this limit in our 
analysis would not qualitatively change our results. 
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explicitly for comparison, but only the LHC2011 data enters our analysis also as part of the 
constraints implemented in HiggsBounds. As explained above, HiggsBounds considers only 
the single channel with the highest expected sensitivity for determining 95% CL (combined) 
exclusion. In the considered region the expected sensitivity of the CMS search [51] happens 
to be slightly higher than the one from ATLAS [50], so that only the CMS limit actually has 
an effect in our analysis. The plot shows also some allowed points with i?^^ above the CMS 
2011 exclusion curve. For these points another channel has a higher expected sensitivity, 
so that the 77 channel has not been selected by HiggsBounds for determining the 95% CL 
limit. 

As one can see in the left column of Fig. [3], for i?^^ in principle a large enhancement, 
roughly up to a factor six, would be possible in the mass range Mh = 114 .. . 130 GeV (and 
an even stronger enhancement for lighter masses). Such large enhancements are now ruled 
out by the LHC searches in the 77 channel. For the points that are allowed by all the 
considered constraints we find that in the region above the LEP limit a suppression of i?^^ 
by more than an order of magnitude is possible. A maximal enhancement of about ~ 3 times 
the SM value, on the other hand, occurs for Mh ~ 124 GeV. This is in interesting agreement 
with the recent data announcing a discovery of a new state compatible with a Higgs boson 
close to ~ 125 GeV, that has been reported both by ATLAS and CMS. This observed excess 
is compatible with a SM Higgs signal, or even better with a signal of a somewhat enhanced 
strength in the two photon channel. Our results show that the MSSM could account for an 
enhanced Higgs signal as compared to the SM case around = 125 GeV with the maximal 
strength that is allowed by the present limits from ATLAS and CMS. The detailed origin 
of this enhancement will be discussed below. Fig. |3] also shows that the possible size of the 
enhancement decreases for larger M^, for Mh = 130 GeV R^^ is confined to values close to 
unity for the allowed points in the parameter space. 

The right column of Fig. [3] shows the corresponding results for the heavy CP-even Higgs. 
For Mh ^ 130 GeV the results for the heavy MSSM Higgs are qualitatively similar to the 
ones for the light CP-even Higgs. In particular, also in this case a slight enhancement over 
the SM rate is possible for Mh ~ 125 GeV for the scan points that are in agreement with 
the collider constraints (as discussed above, the agreement with [g — 2)^ and the observables 
in the flavour sector could be improved by modifying some of the SUSY parameters that do 
not directly influence Higgs phenomenology). The size of the possible enhancement turns 
out to have some sensitivity to the condition that we have imposed on the elements of the 
Z matrix to ensure perturbatively reliable results, see above. Somewhat larger enhancements 
would be possible if the upper limit that we have imposed on the relative size of the two- 
loop corrections were relaxed. Our results for i?;^ demonstrate that the discovery of a new 
boson in the 77 channel at a mass of about 125 GeV that was observed by ATLAS and CMS 
could also be interpreted within the MSSM as arising from the heavier CP-even Higgs boson, 
as discussed in [SU]. Such a scenario would imply that besides a possible signal at about 
125 GeV there would be a lighter Higgs in the spectrum, having significantly suppressed 
couplings to gauge bosons. For Mh ^ 135 GeV we always find i?;^ < 1, in accordance with 
the decoupling properties of the MSSM. 

Before turning to a discussion of the specific mechanisms responsible for possible enhance- 
ments of the di-photon channel in the MSSM as compared to the SM case, we would like to 
make a brief comparison of our MSSM results to those existing in the literature. A previous 
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study of the di-photon channel in the MSSM [68] found a smaller maximum enhancement, 
i?^^ ~ 1.5, and only for a very specific Higgs mass region. We find a larger enhancement, 
primarily due to a larger scan range in /j,. Restricting our scan ranges accordingly, our re- 
sults agree with those of [68]. Similarly, the analysis of [69] uses a more restricted range of 
/X < 1 TeV, and we also find agreement with our results for R^^ in this parameter region. 
In [To], it was claimed that a light Higgs boson with Mh ~ 125 GeV and R'^^ > 1 cannot be 
realised in the MSSM. Clearly, the results from our scan do not corroborate these conclusions 
since we find that such a Higgs boson can indeed be realised in the MSSM, possibly even 
with a rate enhanced compared to the SM. 

The issue of a possible enhancement of i?^^ for a Higgs mass around 125 GeV has also 
been discussed in [31], where in particular the contributions from light staus to BIi{h — > 77) 
and the suppression oi h ^ hh due to Higgs mixing effects have been emphasised. As seen 
above, we find that T{h — )■ 77) can exceed its SM value, which is found to be an effect of the 
stau loop contributions. The most sizeable enhancements observed in i?^^, however, mainly 
arise from a suppression of the total width, which in the SM is dominated by the partial 
decay width into hb. Suppressing the hh channel can therefore yield a significant reduction of 
the total MSSM width. Such a suppression can happen in two different ways. The reduced 
hhb coupling in the MSSM is given at tree-level by 

9hbb sin a ^^^^ 



9HsMbh cos/3 

where a is the mixing angle in the CP-even Higgs sector. In the decoupling limit {Ma ^ 
Mz) the SM is recovered, i.e. (—sin a/ cos/3) — )■ 1. Higher-order contributions from Higgs 
propagator corrections can approximately be included via the introduction of an effective 
mixing angle, corresponding to the replacement a — )■ a^s [71] (in our numerical analysis we 
treat propagator-type corrections of the external Higgs bosons in a more complete way, which 
is based on wave function normalisation factors that form the Z matrix [I?]). A suppression 
of the h ^ hh channel thus occurs for small acs- This feature motivated the definition of 
the "small aefr" benchmark scenario, see [72]. In this scenario, the suppression of r(/i — hb) 
occurs for large tan/3 and very small values of M^, below 200 GeV. 

Genuine corrections to the hhh vertex can lead to another type of suppression. Beyond 
leading order, loop-induced Yukawa couplings of h quarks to the "wrong" Higgs doublet 
are induced. The modified hbb coupling can then be expressed as 

9hbb 1 [ sinacflF , ^^cosOeffA ^^^^ 



9HsMbb 1 + A;, V COS /3 sin /3 

Via the quantity A;, [iniEll terms of C((as tan/3)") and C((ai tan/3)") can be resummed. 
The most relevant contributions are given by 

Ab = ^^^^^ tan/3mg/i J(m?^,m|,m?) + ^^^^^ tan/3 A At /("^|, "^|, l/ip), (33) 



with 



ahlnih/a) -\- achiialc) + hc\Yiic/h) , , 

/ a, 6, c = — ^ . 34 

(a — c)(c — b)[b — a) 
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Figure 4: Results from the MSSM parameter scan on the ratio R^^ in the (MA,tan/3) 
plane. The grey points are excluded by the limits from direct Higgs searches (from 
HiggsBounds 3.6.1, i.e. including LHC2oii)- MSSM Points with an enhancement of i?^^, 
corresponding to > 1, are indicated in black, points with i?^ > 2 are shown in red. 



The dominant higher-order contribution to Aj, are the QCD corrections, given in [75]. Those 
contributions are not included in our analysis, but their effect can be approximated by using 
a scale of rrit for the evaluation of the one- loop expression, Eq. ( 133|) . While the loop-corrected 
hbb coupling, Eq. fl32|) . approaches the tree- level coupling, Eq. fj3T|) . in the decoupling limit 
{Ma ^ Mz), a suppression of g,^f,i is possible for not too large Ma if A;, is numerically sizable 
and positive. For fi > 1 TeV we find enhancements of i?^^ of more than 1.5 for values of Ma 
up to roughly 450 GeV and with moderate to large values of tan/3. Points with i?^^ > 2 
are possible if A;, is relatively large. A;, ^ 0.5. The corresponding effect on i?^^ can be seen 
in Fig. m for the {Ma, tan (3) plane. The points with i?^^ > 1 are indicated in black, and 
the ones with i?^^ > 2 are shown in red. The regions with only grey points are excluded 
by the limits from the Higgs searches at LEP and the LHC2011, where the latter ones arise 
from the searches for MSSM Higgs bosons in the r+r~ channel jlEj. As can be seen from 
the figure, a slight enhancement of i?^^ > 1 can be accommodated over the whole allowed 
region displayed in Fig. |U while scenarios with i?^^ > 2 tend to be closer to the boundary 
of the region allowed by the LHC Higgs limits. 

We now turn to the alternative case where h is light and has suppressed couplings to 
gauge bosons, whereas the heavier CP-even Higgs if is a SM-like Higgs boson. One finds a 
similar enhancement for R^^, which is due to the suppression of gnbiy if 

9Hbb ^ cos a 
9HsMbb cos/3 

is small. Such an enhancement is restricted to the mass region Mh ^ 130 GeV, since for 
higher mass values the coupling of the heavy CP-even Higgs to gauge bosons is suppressed, 
so that the partial width r{H — )■ 77) is smaller than for the SM case, see Fig. [2l Accordingly, 
the scenarios with i?:^ > 1 are only realised in a relatively small parameter region close to 
the exclusion bounds from the Higgs searches, for Ma ^ 150 GeV and intermediate tan/3. 
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Figure 5: Results from the MSSM parameter scan on the ratios R^w hght CP-even 

Higgs bosons h (left column) and -R^^ for the heavy CP-even Higgs boson H (right column). 
The plots in the upper row show R^^y and R^yy as a function of the respective Higgs mass. 
The 95% CL exclusion limits for the WW'^*^ channel using LHC2011 data from ATLAS [52] 
(solid line) and CMS |53] (dashed line) are also shown. The plots in the lower row show the 
correlation to -R^J: and Ri^ . The colour coding is the same as in Fig. [2l 



The scenario in which i?^ is enhanced is complementary to the one giving an enhancement in 
i?^^ (as we checked explicitly). Consequently, a simultaneous enhancement in the di-photon 
channel for both CP-even Higgs bosons is not possible. 

A reduction of the total width, by the suppression of the h,H bb channel, can also 
affect the search for the Higgs boson in other channels. In Fig. [5] we investigate the decays 
h,H ^ WW^*\ The colour coding is as in Fig. [21 but here the green lines show the 95% 
CL exclusion limits for the data presented in 2011 for the channel from ATLAS [52] 

(solid) and CMS [53] (dashed). Here we do not show new results from [1], since this channel 
was updated only by CMS, and no combination of 2011 and 2012 data is available. The 
2012 exclusion curve of CMS is very similar to the one based on 2011 data. Similarly to 
i?^ , also R'^iY can be enhanced in the MSSM. The maximal enhancement of R^^y ~ 2.5 
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At /MsusY At/MsvsY 

Figure 6: Results from the MSSM parameter scan on the ratios i?^^ (left) and i?;^ (right) 
as function of At /Msusy- The colour code is as in Fig. [2J 



that is compatible with the existing bounds is found for ~ 125 GeV. These points, 
above the green lines, are not excluded by the h — )■ H/iy^*) limits, because another channel 
has a higher expected sensitivity for those points (and consequently the channel 
has not been selected by HiggsBounds for determining the 95% CL limit - which would 
not change by the inclusion of the 2012 data shown in [1]). Since the main contribution to 
— )■ 77 comes from diagrams containing VT-boson loops, a strong correlation between these 
two channels is expected. This correlation is confirmed by plots in the lower row of Fig. |5l 
Contrary to T(h — )■ 77), T(h — >■ VTVr^*)) does not exceed its SM value, and therefore the 
maximal enhancement of R^w always smaller than the enhancement of i?^^. This leaves 
some room to have e.g. a SM-like rate, or even a rate that is somewhat suppressed w.r.t. 
the SM, for h — )■ H^IV*-*-* (and h — )■ ZZ^*^) with a simultaneous enhancement in /i — )■ 77. 
Concerning the heavy CP-even Higgs boson, the results are qualitatively similar, but with 
a smaller possible enhancement of the 77 or WW^*^ rate, as we show in the right column of 
Fig. [51 

Finally, in the left (right) plot of Fig. |6]we show the dependence of i?^^ (-R^) as a function 
of At /Msusy, with the colour coding as in Fig. [21 Here it should be kept in mind that values 
of |Af I > a/6 Msusy could potentially lead to charge or colour breaking minima [76]. For the 
light CP-even Higgs, as shown in the left plot, the largest values of i?^^ are indeed found 
around At ~ v^MgusY- However, even applying a stringent cut \At\ < \/6MgusY would 
still leave allowed points with i?^^ ~ 2. For the heavy CP-even Higgs, we see that a stringent 
cut on \At\ would leave nearly all points with i?:^ > 1. 

4.2 CP-even Higgs decays in the NMSSM 

We now turn to the NMSSM and analyse the di-photon decay in this model. For a light 
CP-even Higgs boson with a mass below 114.4 GeVo it was pointed out in [77] that a strong 

As before, we neglect here and in the fohowing, the theory uncertainty of the Higgs boson mass 
evaluation, which for the light Higgs boson should be comparable to or slightly larger than the respective 
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enhancement of the 77 rate — up to a factor of seven over the SM rate — is possible. In 
hght of the recent results, the main focus of our analysis will naturally be on higher Higgs 
masses. For other approaches regarding an interpretation of a possible Higgs signal around 
Mh = 125 GeV in terms of an NMSSM Higgs boson, see [33H35]. 

As before we consider the one-loop induced Higgs decay hi 2 77, but now calculated 
using the framework described in Sect. [31 We perform a scan over the NMSSM parameter 
space and evaluate the partial widths and branching ratios for this mode. The parameter 
ranges used for the scan are given in Tab. [3l The remaining parameters are fixed to the 
benchmark scenario defined above, see Eq. f l25|) . It should be noted that the ranges in 
Tab. [3] are not meant to cover the full NMSSM parameter space. The effects discussed above 
that can cause an enhancement of -R^:^^ in the MSSM can be realised also in the context of 
the NMSSM. In the present analysis we are interested in genuine NMSSM effects, which go 
beyond the MSSM phenomenology. Such genuine NMSSM effects arise in particular from 
the mixing of the Higgs doublet fields with the Higgs singlet. To be specific, we consider 
scenarios that are characterised by large values of Mh± , corresponding to the "SM+singlet" 
limit of the NMSSM discussed above. We furthermore restrict /Xes and tan/3 to relatively 
small values, while our MSSM scan (cf. Tab. [2]) had extended to rather large values of /i and 
tan/3 and had focussed on the region of relatively low values of M^. The parameters are 
chosen such that the mechanisms for enhancing R^!^:^ realised in the MSSM do not play a 
role, putting the emphasis on the genuine NMSSM effects. 



Parameter 


Minimum 


Maximum 




At = Ab = Ar 


-2 400 


2 400 


GeV 




150 


250 


GeV 


Mh± 


500 


1000 


GeV 


tan /3 


2.6 


6 




A 


0.5 


0.7 




K 


0.3 


0.5 




A. 


-100 


-5 


GeV 



Table 3: Parameter ranges used for the CP-even Higgs decay scan in the NMSSM. 

The results for hi^2 77 are shown in Fig. [7] [h^ is always heavy and plays no role in 
our analysis). The colour coding is the same as in Fig. [21 i.e. all displayed points satisfy 
the "theoretical" constraints of Sect. 13.4] (grey). Points which in addition fulfil the direct 
Higgs exclusion limits from colliders (from HiggsBounds 3.6.1, i.e. including LHC2011) are 
drawn in blue, and points which satisfy all the constraints, in particular also those from Aa^ 
and flavour physics, are shown in black. The red curve shows the corresponding SM result, 
obtained by setting Mh^^ = mh^^^^, respectively. 

We choose to study m/u in a range from 80 GeV up to its maximum around 135 GeV. Al- 
lowed points with < 80 GeV are also found in the scan, but the large singlet component 
of these very light Higgs bosons gives rise to a quite different phenomenology, which we do 
not investigate in detail here. For masses close to 140 GeV, the number of allowed points is 
seen to decrease, which illustrates that only quite specific choices of the input parameters give 

uncertainty of the MSSM, i.e. roughly at the level of 2 - 3 GeV [211. 
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Mh, (GeV) Mh, (GeV) 

Figure 7: Results from the NMSSM parameter scan (see text) for the partial widths T{hi — >■ 
77) and the corresponding branching ratios of hi (left) and h2 (right). All points in the 
figure fulfil the "theoretical" constraints defined in Sect. 13.41 In addition, the blue (dark) 
points satisfy direct Higgs search limits from colliders (from HiggsBounds 3.6.1, i.e. including 
LHC2011), while the black points are in agreement with all theoretical and experimental 
constraints. The solid (red) curve shows the respective quantities evaluated in the SM. 

close to the maximum. This, as well as other features with local under- (over-) density 
of points in certain regions, can simply be viewed as sampling artefacts, i.e. the point density 
has no physical meaning. For h2 we study the mass interval 120 GeV < < 170 GeV, 
which means there is an overlap with the region considered for m/^ . To go even higher in rrih^ 
is not particularly interesting for our purposes, since when the two-body decay /i2 WW^*^ 
is open the loop-induced h2 — )■ 77 decay becomes suppressed (as is also clearly visible in the 
figure). 

Fig. Hshows that T{hi — 77) is always smaller than (or at most equal to) its SM value for 
the points in our scan. This means in particular that our scan, for which we have fixed the 
slepton masses to large values (see Eq. fl2Sl) ). does not contain points with light staus (the 
contribution of light staus was discussed in the MSSM context above). For ^ 140 GeV, 
the partial width does not reach the full SM value, which shows that this mass region is 
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Figure 8: Results from the NMSSM parameter scan on the ratio B!^ty for the two hghtest 
Higgs bosons hi (left column) and /i2 (right column). The plots are displayed both on a 
logarithmic scale (upper row) and on a linear scale (lower row). The colour coding for the 
scan points is the same as in Fig. [71 The green lines show exclusion limits on this channel 
at 95% CL from LHC2011 data from ATLAS [50] (solid) and from CMS [5l] (dashed). The 
red lines are the new limits from ATLAS (solid) and CMS (dashed) taken from |T]. 



not accessible for a fully SM-like Taking into account the collider constraints, we also 
see that, as in the MSSM, a SM-like Higgs boson with m^^ < 114.4 GeV is excluded as a 
consequence of the LEP limits. Despite the smaller NMSSM width for T{hi — )■ 77) compared 
to the SM, Fig. [7] shows that an enhancement of the branching ratio with up to an order of 
magnitude over the SM is possible. The results are similar for hi and h2 in the overlapping 
mass region. 

As in the case of the MSSM we now analyse Bl^i^. The total widths appearing in Eq. (l30l) 
are calculated in an approximate way according to 

+ T{h ^ WW^*^) + T{h -> 77) + r(/i. ^ gg). (36) 

where T^h^^. denotes the one loop self energy of h^. The inclusion of the off-shell decays, as 



^tot{hi] 



Tm 
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Figure 9: Results from the NMSSM parameter scan on R^i^ and the elements U^j of the 
CP-even Higgs mixing matrix for hi (top) and /i2 (bottom). The colour coding is the same 
as in Fig. [71 



well as the loop-induced processes, in the total width is essential for a realistic prediction. 

The results for R^l^ and R^l^ from the scan over the NMSSM parameter space are shown 
in Fig. IHl As before, we show the plots both on a logarithmic and a linear scale, and 
include the latest results on this channel taken from [1]. Looking first at hi, the figure 
shows that a sizable enhancement over the SM rate is possible over the whole mass range 
from TTth^ = 80 GeV to rrihi = 130 GeV. For the range of Higgs masses below the SM 
limit, m/i^ < 114.4 GeV, points with a significant enhancement i?^;^ ^ 7 are observed, in 
accordance with the results of [77] (see also [5B]). 

Turning to /12, the results for i?^^ are similar to those for R!]jl^ in the common mass 
range; the observed maximal enhancement is i?!^^ ^ 2 for m/^j in the range from 120 GeV 
to 125 GeV. A smaller enhancement over the SM is possible for all < 145 GeV. As 
approaches 160 GeV, where the on-shell decay /i2 — > VTty^*) opens, the rate drops to 

Comparing the results for Ri!^i^ to the limits most recent limits from ATLAS (solid red) 
and CMS (dashed red) [1] (where the LHC2011 limits are shown in green), it is clearly visible 
that the NMSSM (similarly to the MSSM) can produce points with a large suppression of 
Ry^. Concerning the more interesting case i?;!;; > 1, we see that the NMSSM can produce an 
enhancement compatible with an excess over the SM rate for Higgs production in the mass 
region around 125 GeV. The results of Fig. [8] show that such observed excess over the SM 
rate for Higgs production in the 77 channel is well compatible with both hi or /i2 production 
in the NMSSM. 

In order to identify the conditions under which a significant NMSSM enhancement of 
R^* is possible — as explained above, this is a genuine NMSSM effect that goes beyond the 
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mechanisms discussed above for the MSSM — we show in Fig. [9] the dependence of R^i^ on 
the composition of hi as defined in Eq. f[T^ . Fig. IH] shows that an important requirement 
for R'^i^ > 1 is that ^ 0, which means that the corresponding Higgs mass eigenstate 
lacks a Hd component. In the hmit — ?■ (corresponding to ttofr — in the MSSM) the 
otherwise dominant decay channels hi — hb and hi — ?■ rr vanish (in the effective coupling 
approximation), thereby increasing BR(/ij — 77). Since is unitary, the general sum rule 

El^fl' = l (37) 
3 

implies that points with = must have |?7jf p + |f/jf p = 1. From Fig. [9] it can be seen 
that a configuration that maximises R!^i^ would be 

Unlike the case of Higgs doublet mixing resulting in a small in the MSSM (which requires 
a low value for Ma and a high /i), in the NMSSM this feature is caused by a sizable singlet 
component of hi. The observed R'^iy enhancement is therefore a genuine feature of the 
NMSSM which is still present even in the SM+singlet limit. On the other hand, in the MSSM 
limit (where p = 0) points from our scan show only very small -R^j^ enhancements. This 
is a consequence of the fact that we have restricted our scan in the NMSSM to large Mu±, 
large slepton and squark masses, as well as to relatively small values of tan (3 and /lefr, which 
corresponds to a parameter region in the MSSM that is complementarity to the one used for 
our MSSM scan. 

As for the MSSM case, also in the NMSSM the decay hi — )■ 77 is usually dominated by 
contributions from loops containing W bosons, and we expect a corresponding correlation 
of i?^^ with the (off-shell) decays hi WW^*^ and hi ZZ^*\ This is studied in Fig. [IHl 
where we give results for the tree-level decays hi^2 WW'^*'^ using the same colour coding 
as above. As expected, a sizable enhancement is possible for Ry^ry/, in particular for hi, 
and a strong positive correlation between R!Iyw ^^"^ -^77 visible. The possibility of a 
simultaneous enhancement of these two modes can again be understood as an effect of the 
large suppression of the main fermionic coupling /ijfefe, which leads to an enhancement of the 
respective branching ratios. As in the case of the MSSM, despite the positive correlation 
between R^i^ and R^i^ it is nevertheless possible to have a both a slight enhancement of 

R^i^ and a slight suppression of Ry/w- 

In Fig. do] we also display the ATLAS [52] and CMS [53] limits on the H WW channel 
based on LHC2oii- For the same reasons as in the MSSM case we do not include the new 
results from [1]. The LHC2011 limits exclude a SM-like Higgs for Mh ^ 130 GeV (which 
would not change by the inclusion of the data from pj). The points excluded by this mode 
can be seen as a grey strip for rrih^ ^ 130 GeV in this figure and also in Fig. |8l For the few 
points above the green lines which are not excluded by — )■ WW, another channel has a 
higher expected sensitivity. For m/j. < 130 GeV the existing limits on hi — )■ WW do not pose 
any constraint on the NMSSM, but an improved exclusion of the WW mode at lower Higgs 
masses could have severe implications for the scenarios we consider. The resulting limits on 
are competing with those from direct searches in this channel, but appear somewhat 
relaxed in view of the excess over the SM prediction in the 77 channel as reported in [1]. 
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Figure 10: Results from the NMSSM parameter scan on the ratios R^J^yy ii = 1,2) for the 
NMSSM Higgs bosons hi (left column) and h2 (right column). The plots in the upper row 
show R'^iy as a function of the respective Higgs mass. The 95% CL exclusion limits for 
the WW^*^ channel from LHC2011 data from ATLAS [52] (sohd line) and CMS [53] (dashed 
line) are also shown. The plots in the lower row show the correlation of Ry/y^ and R^ty. The 
colour coding is the same as in Fig. [71 

The combination of those two channels offers an excellent opportunity for cross-checks and 
model testing. 

4.3 CP-odd Higgs decays in the NMSSM and the MSSM 

For completeness, we have also investigated decays of CP-odd Higgs bosons in the MSSM 
and in the NMSSM. In the following we denote a generic CP-odd Higgs boson as (pA- As a 
consequence of gauge invariance, the '^(PaZ vertex vanishes in both models at the tree-level. 
The vertices VV(f)A iY = 7, 9, Z, W) vanish by virtue of CP-invariance [78]. Accordingly, all 
these vertices only appear at the one-loop level. 

The MSSM contribution to A — )■ 77, gg has been evaluated by many groups (see [18] for 
a comprehensive list of references). The MSSM prediction for A — )■ WW was first studied 
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in Refs. [79l|80], however only heavy SM fermions were included. In Ref. [80], chargino and 
neutralino contributions to A — j- WW were included together with the SM fermions. In the 
MSSM the BR's of A ^ 77 and A WW are of 0(10-^) and C(10-3), respectively, for 
small tan /3 and decrease for large tan /3. Consequently, for a production cross section similar 
to a SM Higgs boson with the same mass, R^^ can be order of (9(10"^). The prospects for 
the observation of those modes at the LHC are not very encouraging. 

In order to assess the NMSSM potential for the observation of the CP-odd Higgs bosons 
in these loop-induced decay channels, a random scan over the six-dimensional parameter 
space of the Higgs sector has been performed (see Tab. |3]). In this scan the lightest CP- 
odd Higgs boson has a mass ranging from about 80 to 250 GeV. Due to the mixing of 
doublet and singlet states, Eq. f|T3|) . the NMSSM CP-odd Higgs bosons can have couplings 
to fermions that are significantly smaller than in the MSSM. Contrary to our analysis in the 
previous section where we compared to the SM, we now compare the NMSSM results for the 
light CP-odd Higgs boson ai to the results for the CP-odd Higgs boson in the MSSM (with 
TBai = Ma)- Both in the MSSM and the NMSSM, we compute the total decay widths of 
CP-odd Higgs according to 



-Im 



J2 ^i^A ^ VV). {31 

V=-y,g,Z,W 



In order to compare the NMSSM to the MSSM we define the ratio of di-photon cross-section 
normalised to the MSSM rate as follows 



pa, ^ Hag ^ flj) X BR(a, ^ 77)]nmssm /oq^ 
- K^?^ ^ A) X BR(A ^ 77)]mssm ■ ^ ^ 



In the case where A, k — )■ (with k/A held constant), which corresponds to the MSSM 
limit (see above), the interactions of the singlet field with the MSSM Higgs doublets vanish, 
i.e. the singlet field decouples. The decoupled singlet field itself has no gauge couplings. In 
this MSSM limit, one finds that one of the NMSSM CP-odd Higgs bosons, au (which does not 
necessarily coincide with ai) behaves like the MSSM A boson. The MSSM phenomenology is 
fully retrieved in the Higgs sector, the total decay width of the CP-odd Higgs in the NMSSM 
approaches the MSSM value, and we find ~ 1. However, even in this limit, in some 
cases the distinction of the two models might still be possible via the neutralino sector [81] . 

Far from the MSSM limit, the coupling a^ff, which is proportional to the doublet-singlet 
mixing, gets suppressed and the NMSSM phenomenology can differ from the MSSM. In our 
scan we find an enhancement of the BR(ai — )■ 77) of up to 100 over the MSSM. However, 
also the oi production rate gets suppressed in this case, leading to (at most) i?"^ ~ 2. 
Consequently, as for the MSSM, also ai is expected to remain unobservable in this channel 
at the LHC (and we do not show the plots for this channel). 
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5 Conclusions 



The discovery of a new particle with properties that are compatible with a Higgs boson has 
started a new era, with the goals to reveal whether the new particle is indeed a Higgs boson 
and whether its detailed phenomenology is consistent with the SM predictions or whether 
there is evidence for physics beyond the SM. The observations reported both by ATLAS 
and CMS that the new state has a signal strength in the 77 channel that is considerably 
higher than the one expected in the SM may be a first hint into this direction, although the 
statistical significance of this deviation from the SM prediction is not sufficien at present to 
draw a definite conclusion. In view of the possibility of an enhanced rate in the 77 channel, 
it is nevertheless interesting to investigate the predictions of possible alternatives to the SM 
in this context. We have carried out such an analysis in the present paper, comparing the 
predictions for Higgs boson production in gluon fusion, the main production channel at the 
LHC, and its subsequent decay into two photons in the SM, the MSSM, and the NMSSM. 
We have furthermore analysed the WW* channel, which is strongly correlated with the 77 
mode. 

While for the predictions in the SM and the MSSM we have used the well-known code 
FeynHiggs, for the predictions in the NMSSM we have developed a new framework consisting 
in particular of an appropriate model file for the program FeynArts and the input on masses, 
mixing angles, etc. needed for the numerical evaluation. In this set-up we can make use of 
the highly automated programs FormCalc and LoopTools for the evaluation of the relevant 
observables. In the present implementation NMSSMTools is used for higher-order corrections 
to the masses of the external Higgs bosons. Numerous tests have been performed to verify 
the implementation; among other things we evaluated more than 150 loop-induced processes 
in the NMSSM and checked the results for UV-finiteness. 

We have presented results for Higgs-boson production in gluon fusion and its decays into 
77 and WW* within the MSSM and the NMSSM, normalised to the SM prediction. We 
have analysed in detail possible mechanisms for the enhancement (but also the suppression) 
of those channels in both models. In this context we have investigated in particular whether 
an enhancement of the 77 rate for a Higgs mass of about 125 GeV is compatible with ex- 
isting limits on the parameter space arising from theoretical constraints as well as from the 
limits from direct searches for supersymmetric particles, from the Higgs searches at LEP, 
the Tevatron and the LHC (where we have incorporated the limits based on the data pre- 
sented in 2011), from electroweak precision observables and from flavour physics. Performing 
parameter scans in both models, we have then confronted the points passing all the above 
constraints with the latest results of the Higgs searches in the 77 channel that have just been 
announced by ATLAS and CMS. We have found that an enhanced rate of Higgs production 
and decay to two photons can easily be realised in both models, the MSSM as well as in the 
NMSSM. This holds not only for the lightest CP-even Higgs boson in the models, but also 
for the second lightest CP-even Higgs boson in both the MSSM and the NMSSM. In this 
latter interpretation in both models the lightest CP-even Higgs boson possesses a strongly 
suppressed coupling to gauge bosons and escapes all existing direct searches. 

Within the MSSM we have analysed the mechanisms that can lead to such a enhanced 
77 rate in comparison to the SM prediction. Besides the presence of light scalar taus, 
in particular a suppression of the bb decay mode results in an enhanced 77 rate. This 
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suppression can either be caused by Higgs-boson propagator corrections entering the effective 
mixing angle, or by the so-called A;, corrections. 

Within the NMSSM the above mentioned mechanisms can also naturally be realised, and 
we focused on additional mechanisms that are genuine for the NMSSM. We found that in 
particular the doublet-singlet mixing can result in a substantial suppression of the bb mode, 
resulting again in the desired enhancement in the 77 rate with respect to the SM prediction. 

Finally we have briefly analysed the decays of CP-odd Higgs bosons to two photons. 
While an enhancement with respect to the corresponding rate in the MSSM is possible in 
the NMSSM, a signal such as reported by ATLAS and CMS [Ij cannot be accounted for by 
a CP-odd Higgs. 
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